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Abstract 

The dynamics of entanglement between two photons with one of them passing through noisy 
quantum channels is characterized. It is described by a simple factorization law which was first 
theoretically proposed by Konrad et al. [Nature Phys., 4, 99 (2008)]. Quantum state tomogra- 
phy process is employed to reconstruct the reduced density matrixes of the final states and the 
corresponding concurrences are calculated. Good fittings between experimental results and the- 
oretical predictions are found, which imply the validity of the general factorization law in the 
characterization of entanglement dynamics. 
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Entanglement is the crucial resource for quantum communication and computation 
However, due to unavoidable couplings with the environment, it is readily destroyed when 
entangled particles transmitted through noisy quantum channels j^J. Characterization the 
dynamics of entanglement becomes an important task for those entanglement-based quan- 
tum information processing protocols and it is also crucial for understanding the distinct 
)roperties of entanglement dynamics, such as the phenomenon of entanglement sudden death 

In order to quantify the entanglement, we need a convenient quantity to measure it. The 
Wootter's concurrence |8[] is a good choice for two qubits, which is given by 

C = max{0,r}, (1) 

where T = a/Ai — — — V^4 and the Aj are the eigenvalues in decreasing order 
of the matrix p(a y ® o y )p*(o y ® a y ). o y is the second Pauli matrix and p* is the complex 
conjugate of p in the canonical basis {|00), |01), |10), |11)}. The concurrence ranges from 
(the separated state) to 1 (the maximally entangled state). 

Usually, the tomography process 9| which requires measurement on a complete set of 
observable quantities is used to reconstruct the final density matrix p and the entanglement 



of the state is always nonlinear dependent on it 



10], just as we have shown above. Therefore, 



the knowledge of the dynamics of entanglement can only deduce from the time evolution of 
the state O, 3, [ill, 12, Q, Q. The general result on entanglement dynamics is still a 
great challenge. 

Recently, the work of Konrad et al. gives a general factorization law on the entanglement 
dynamics of the biparticle system under the action of an arbitrary channel on one of the 
component [isj]. The main result they got is the formulation: 

C[(l ® $)| X )(x|] = C[(l ® $)|0+><0 + |]C(| X », (2) 

where \x) is the initial input pure state and \<p + ) is the maximally entangled state. We 
use Clp to represent the left term of the equality and Crp to represent the right term, 
for simplicity. It is clearly shown that the entanglement dynamics in the one-sided noisy 
channel only relate to the dynamics of the maximally entangled state in the channel with a 
factor of the concurrence of the initial state 

Due to the convexity property of concurrence , this factorization law can be generalized 



to the mixed initial state po [li 

C[(l ® $) Po ] < C[(l ® $)|0 + )(0 + |]C( Po ). (3) 

We also use and (Tra/ to represent the left and the right terms of this inequality, 

respectively. This inequality gives the upper bound of the evolved entanglement. 

Here, using an all optical setup, we experimentally characterize the two-photon entangle- 
ment dynamics with one of the photons passing through the phase damping channel or the 
amplitude decay channel. The experimental results show that the time evolution of entan- 
glement is fully characterized by the channel's action on the maximally entangled state and 
we verify the factorization laws. 

Our experimental setup is shown in Figure 1. Ultraviolet (UV) pulses with wavelength 
centered at 390 nm which are frequency doubled from a Ti:sappire laser with a pulse width 
of 130 fs and a repetition rate of 76 MHz pump into a two geometry type I beta-barium- 
borate (BBO) crystals to generate entangled two photons [16j. The polarization of these 
UV pulses is set by the half- wave plate HWPO which can be used to prepared different pure 
input states \x) = a\HH) a ^ + (3\VV) a ,b (a and (3 are the relative amplitudes which are set 
to be real for simplicity and they obey the relationship of a 2 + (3 2 = 1). Quartz plates (CP) 
in both emitted modes a and b are used to compensate the temporal difference between 



horizontal and vertical polarization components in these two nonlinear crystals [171 ]. 

A Fabry-Perot (F-P) cavity followed by quartz plates Q2 in the dotted pane M is used to 
simulate the phase damping channel, which has been used in our previous experiment 
When the photon passes through Q2 with thickness L, the relative phase between H- and V- 
polarization photon can be calculated as <p = LAnoo/c, where c is the velocity of the photon 
in the vacuum. An = n Q — n e is the difference between the indexes of refraction of ordinary 
light (n ) and extraordinary light (n e ), which can be treated as a constant of 0.01 for the 
small frequency distribution. The amplitude decay channel shown also in the dotted pane 
M is simulated by a sets of glass slabs (gs) which are tilted to the Brewster angles (about 
57°) [18(]. After passing through the slaps, the photon with vertical polarization has some 
probability to reflect and the photon with horizontal polarization will transmit completely. 
Therefore, the corresponding quantum map [3] can be written as 

\H) s \0)e-^\H)s\0)e, 



\V) s \0)e — > VT^~e\V) s \0) E + y/i\H) s \l) E . (4) 
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FIG. 1: (Color online). The experimental setup to investigate the entanglement dynamics in 
different noisy channels. Entangled polarization photon pairs generated from the adjacent nonlinear 
crystals (BBO) emit into modes a and b. The half-wave plate (HWPO) is used to change the 
polarization of the pump light. A half-wave plate (HWP1) and quartz plates (Ql) are inserted 
into mode a depending on different cases. The dotted pane M represents the noisy channels. The 
Fabry-Perot (F-P) cavity followed by quartz plates (Q2) simulates the phase damping channel, 
while the sets of glass slabs (gs) which are tilted to the Brewster angles represent the amplitude 
decay channels (the dashed lines represent the reflected photons). Quarter- wave plates (QWP), 
half-wave plates (HWP) and polarization beam splitters (PBS) in both arms are used to set the 
detecting polarization bases for the state reconstruction. Both photons are finally detected by 
silicon avalanche photodiodes (Dl and D2) equipped with 3 nm interference filters to give coincident 
counts. 

where |0) e and |1) e represent the propagation paths of the photon and e represents the total 
reflectivity of the glass slaps. We can change the polarization of the reflected photon into 
\H) with half- wave plates and detect them with single photon detectors (not shown in fig. 

ED. 

A half-wave plate (HWP1) and quartz plates (Ql) are inserted into mode a depending 
on different cases to introduce different decoherence effect on the maximally entangled state 
\(f) + ) = -^{\HH) a £ + \ VV) a ,b) to prepare different mixed input states. The optical axis of 
HWP1 is set to be 22.5° operating as a Hadamard gate and Ql with the optical axis set to 



4 



be horizontal can induce relative phase in the basis H/V. 

The density matrixes of the final states are reconstructed by the tomography process 
9|]. Experimentally, quarter- wave plates (QWP), half- wave plates (HWP) and polarization 
beam splitters (PBS) in both arms are used to set the standard 16 polarization analysis 
measurement bases [9(. Both photons are detected by silicon avalanche photodiodes Dl and 
D2 to give coincident counts. Narrow band interference filters (IF) with a full width at half 
maximum (FWHM) of 3 nm in front of detectors are used to reduce the background and 
define the bandwidth of the photons. 
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FIG. 2: (Color online). Experimental results for the entanglement dynamics of different pure initial 
states in the phase-damping channel. (a)a 2 = |. (b) a 2 = |. (c) a 2 = j^. Red dots and blue 
squares represent Clp and Crp, respectively. The solid lines and dashed lines representing the 
theoretical predictions of Clp and Crp completely overlap. Ao=780 nm. 

We first demonstrate the entanglement dynamics of different pure input states in the 
phase-damping channel, as shown in fig. 2. Three different pure input states \\) with 
a 2 = | (fig. 2(a)), a 2 = | (fig. 2(b)), a 2 = (fig. 2(c)) are considered, respectively. 
The x axis represents total thickness of Q2 which is given by its retardation and Ao = 780 
nm. Experimental results of Clp and Crp (the product of C[(l £g> $)|c/> + ) (</> + |] and C(|x))) 
are represented by red dots and blue squares, respectively. We can find that both of them 
overlap with each other very well, which infers that Clp = Cpp. The theoretical predictions 
Clp (solid lines) and Crp (dashed lines) are calculated from equation ([TJ) with the fitting 
parameters the same as the case we demonstrate the phenomenon of entanglement sudden 
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FIG. 3: (Color online). Experimental results for different mixed input states, (a) the case with 
inserting only Ql=117Ao into mode a. (b) the case with inserting both the HWP1 and Ql=117Ao 
into mode a. Red dots is the experimental results of Tlm and blue squares represent Crm- The 
solid lines and dashed lines are the theoretical predictions of Clm and Crm, respectively (they 
completely overlap in the case of (a)). At the area of Tlm < where the theoretical prediction is 
represented by dotted lines (b), the concurrence is set to according to equation (1). 

death [?J and they completely overlap with each other. It is also seen from fig. 2 that 
experimental results agree well with theoretical predictions. In our experiment, the pumping 
power is about 200 mW and the integration time is 6 minutes. Error bars are mainly due 
to the counting statistics and the uncertainties in aligning the wave plates 9]. 

When it cames to the case with mixed input states, the dynamics of the maximally 
entangled state in the channel gives the upper bound of the entanglement dynamics according 
to inequality fl3]). Fig. [3^a) shows the case with the mixed input state prepared by only 
inserting Ql=117Ao into mode a to dephase the maximally entangled state \<fi + ) {a 2 = 
|). We can see that the experimental results agree well with the theoretical prediction 
and Clm = Crm- We further demonstrate the other case in which we insert both the 
HWP1 and Ql = 117Ao- Fig. [3]^b) shows our experimental results, where the phenomenon of 
entanglement sudden death {Clm) occurs [7|. At the area of Tlm < 0, the concurrence is 
set to corresponding to equation (1). It can be seen that Clm (red dots) is less than Crm 
(blue squares) during the evolution, which is consistent with the theoretical predictions of 
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Clm (solid line) and Crm (dashed line). 




FIG. 4: (Color online). Entanglement dynamics of different pure states in the amplitude decay 
channel, (a) a 2 = |. (b) a 2 = |. (c) a 2 = j^. Red dots and blue squares represent experimental 
results of Clp and Crp, respectively. The corresponding theoretical predictions of Clp (solid 
lines) and Crp (dashed lines) completely overlap. 

Now, we consider the entanglement dynamics in the amplitude decay channel, where the 
F-P cavity and quartz plates Q2 in the dotted pane M are replaced by the set of glass slaps 
titled to the Brewster angles. By changing the number of the glass slaps, we can control the 
reflectivity of this channel. In order not to disturb the transmitted photons, two relatively 
placed slaps are added or removed at the same time. The reflectivity of each two slaps 
we measured is about 0.46. The counts of reflected photons is deduced from the difference 
between the transmitted vertical polarization photons without glass slaps and with glass 
slaps. The coincident counts of the 16 measurement bases to reconstruct the final photon- 
state are consisted by the corresponding transmitted and reflected parts. The transmitted 
counts of each measurement basis are detected by the tomography measurement setup in fig. 
[TJ while the corresponding reflected counts is calculated as the product of the total counts 
of the reflected photons and the theoretical detecting probability of each measurement basis 
in the reflected photon-state. 

Fig. H] shows the case of the amplitude decay channel with different pure input states |%). 
In this experiment, the reflected photon-state is \VH) a ^. Fig. BJa), (b) and (c) correspond 
to the initial state with a 2 — |, | and ^, respectively. We can see that the experimental 
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results of Clp (red dots) are equal to Crp (blue squares). The theoretical predictions of Clp 
(solid lines) and Crp (dashed lines) are both given by 2a/3\/l — e. The difference between 
the experimental results and the corresponding theoretical predictions comes mainly from 
disturbance of the preparation of the initial state and the absorbtion of the slabs. 
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FIG. 5: (Color online). Entanglement dynamics of the amplitude decay channel for mixed input 
states, (a) the case with inserting only Ql=117Ao into mode o. (b) the case with inserting the 
HWP1 and Ql=117Ao into mode a at the same time. Experimental results are represented by red 
dots (Clm) and blue squares (Crm)- Solid lines and dashed lines are the corresponding theoretical 
fittings of Clm and Crm (they completely overlap in the case of (a)). 

Fig. [5] shows the experimental results of the amplitude decay channel for mixed input 
states. The initial mixed state in fig. Ufa) is the same as that in the case of fig. 3(a), which 
is prepared by inserting only Q1=117A into mode a with a 2 = ~. We find that Clm (red 
dots) equal to Crm (blue squares). The theoretical predictions of Clm (solid lines) and 
Crm (dashed lines) completely overlap and are consistent with the experimental results. 
Therefore, the entanglement dynamics of the inequality ([3]) reaches its upper bound. In 
another case where the initial mixed state is prepared by inserting the HWP1 and Ql=117Ao 
at the same time, we find that Clm < Crm and inequality (3) holds, as shown in fig. [5](b). In 
this experiment, the reflected photon-state is given by p r = 1/2(\HH)(HH\ — k*\HH)(VH\ — 
K a \VH) (HH\ + \VH) (VH\), where K a is the decoherence parameter in mode a. Experimental 
results agree with the corresponding theoretical prediction employing equation flTJ). 
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In conclusion, we have characterized the entanglement dynamics in quantum channe 



with the maximally entangled state and verified the entanglement factorization law 



15| 



The results provide a novel method to describe entanglement dynamics in noisy channe 



and would have great importance on the construction of the complex quantum network 



20] 



Note: After we finished this manuscri pt, we find that similar experimental results have 
been published online by Farias et al. [21]. They have characterized the entanglement 
dynamics in the amplitude decay channel, which is simulated by a modified Sagnac interfer- 
ometer. In our experiment, we demonstrate the characterization both in the phase damping 
channel and the amplitude decay channel, which are simulated by the F-P cavity followed 
by quartz plates and the glass slaps with the axes titled to the Brewster angles, respectively. 
What is more, the concurrence of the evolved maximally entangled state (C'[(1®$)[^")(0 + |]) 
is calculated from the experimental form of the single channel $ in their paper 21(, while 
we measure all the quantities in equality ([2]) and inequality ([3]), and then compare the left 
terms and the right terms to verify the idea of characterizing the entanglement dynamics in 
noisy channels with the maximally entangled state. 
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